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Abstract 

We discuss the properties of interacting electrons on a finite chain with open 
boundary conditions. We extend the Haldane Luttinger liquid description to 

o ■ 

■rj" ' these systems and study how the presence of the boundaries modifies various 

o 

■ correlation functions. In view of possible experimental applications to quan- 

"j^ ■ tum wires, we analyse how tunneling measurements can reveal the underlying 

Luttinger liquid properties. The two terminal conductance is calculated. We 

^ . also point out possible applications to quasi one dimensional materials and 

O 



study the effects of magnetic impurities. 

•i-H 

X 

: I. INTRODUCTION 



Physical properties of one dimensional (ID) metals are by now well understood theoret- 
ically. Unlike higher dimensional metal, where interaction slightly modifies the free Fermi 
gas behavior (Landau- Fermi liquid picture), in ID metals the electron-electron interaction 
plays a fundamental role and strongly affects their physical properties as compared to the 
ideal Fermi gas. The most interesting feature is the absence of well defined single particle 
excitations. The only stable low energy excitations turn out to be collective charge and spin 
density fluctuations (zero sound modes). Charge and spin sounds are dynamically indepen- 
dent, which, together with the absence of the continuum of electron-hole excitations, gives 
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rise to the so-called spin-charge separation. The interaction also modifies the asymptotic 
behavior of all the correlation functions: at large distances (times) they are shown to decay 
as a power law with interaction dependent exponents (Luttinger liquid behavior, for a review 
see Ref.0). 

Before the recent achievements of the submicron-size technology in the fabrication of ID 
quantum wires, ID electron gases have been studied either on their own rights (as interesting 
mathematical objects) or in view of applications to quasi ID materials!. For these purposes, 
it was sufficient to investigate infinite systems (or to impose periodic boundary conditions, 
which are relatively easy to treat). Semiconductor devices with ID confining potentials 
(quantum wires) represent a new and promising realizations of ID electron systems which 
provide an alternative way to study experimentally the Luttinger liquid properties. For 
instance, the collective nature of the low energy excitations have been successfully probed 
by inelastic light scattering experiment sH. Other experiments particularly suited for quantum 
wires are those measuring transport properties, and an interesting question is whether these 
measurements are able to reveal the Luttinger liquid character of ID systems. As regards 
to the measurements of bulk properties in a clean system (e.g. optical conductivity) the 
answer would be no. In fact, although many properties are anomalous, the (bulk) transport 
properties of a clean ID metal are expected to be qualitatively similar to those of an ordinary 
3D metal. As regards to the surface measurements or, more generally, to the response to 
any local probe, it is claimed that the answer is positive. This result has been reached by 
extending the analysis of idealized ID chains to physical quantum wires which are finite 
systems with open boundaries or contacts of ID electron gases with normal 3D metals. 

The anomalous response of an interacting ID Fermi gas to a local probe was firstly 
recognized by Kane and Fisher^. Subsequently many different experimental situations have 
been proposed and analysed in the framework of the Luttinger liquid theory of ID metal^. 
Several effects due to open boundary conditions have been studied by means of conformal 
field theoryi. All these studies suggest that the interaction affects (in a non trivial way) the 
behavior of the Fermi gas close to the boundaries. Consequently transport measurements 
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probing the edge properties should in principle reveal the Luttinger liquid behavior. 

In this paper we study the properties of a ID chain of interacting electrons with open 
boundary conditions from a more traditional point of view not assuming conformal invari- 
ance. Some of the results will be new, some not; our intention is to provide a simple but 
powerful tool to tackle different problems that arise while studying finite ID systems. 

The layout of the paper is as follows. In Section II we develop the bosonization tech- 
nique appropriate for open boundary conditions, generalizing Haldane's approach (which 
was originally deviced for periodic systems). The bosonization for non interacting electrons 
is discussed first. Then we study the interacting case. In Section III several correlation func- 
tions are calculated and the Priedel oscillation caused by the boundary is discussed. We also 
give an alternative derivation, based on Ward's identity method, of correlation functions and 
outline some applications of the open boundary analysis to slightly doped quasi- ID electron 
systems. Section IV is devoted to the analysis of transport phenomena in quantum wires. 
The two terminal conductance is calculated and some other tunneling processes at the edge 
of the wire are discussed. In Section V we study a magnetic impurity coupled to a ID elec- 
tron gas in the presence of a potential scattering. There are four appendices. In Appendix 
A some formulae used in computation of correlation functions are given. The remaining 
appendices are devoted to extensions of our approach: an X-ray edge type problem related 
to the scattering potential sign changing operator is investigated in Appendix B, effects of 
long-range electron-electron interaction are examined in Appendix C, and gapped phases of 
finite systems are studied in Appendix D. 

II. OPEN BOUNDARY BOSONIZATION 

This section is intended to give a general method of bosonizing an interacting Fermi- 
system in the case of open boundaries (i.e. of a finite system of the length L). 

The bosonization method has a long history. The equivalence between the excitation 
spectra of interacting fermions and free bosons in one dimension (ID) was established by 
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Mattis and Liebi in their solution of the Luttinger mo dei. The bosons were identified as 
particle-hole excitations over the Fermi sea, and their dynamics turned out to be that one 
predicted by random-phase approximation. However, the full power of the bosonization 
method became clear later, when the representation for the electron creation (annihilation) 
operator in terms of free bosons was discovered&^iil. This provided a powerful tool to calcu- 
late fermionic correlation functions in terms of free boson correlation functions. (Actually, 
this representation is very close in spirit to the famous Jordan and Wigner representation.) 
The finite size effects (assuming periodic boundary conditions) were studied by Haldandi. 
We start with the Hamiltonian: 

H = Hq + Hint , 

where the first term represents the kinetic energy, 

Ho=Y. r dx^l{x)e{-td,)^P,{x) , (1) 
and the second one describes the electron-electron interaction. 

Hint = \ll j dxdy^lj\{x)ijl{y)Uss'{x - y)^ ^,{y)^lj ^{x) , (2) 

ss' 

e{k) is the dispertion law of the ID band, and ips{x) is the spin s electron annihilation 
operator subject to the open boundary conditions: 

^,(0) = ^s{L) = (3) 



A. Non-interacting electrons 

We consider the non-interacting case first. This situation has already been discussed in 
the literature (see e.g. Ref. [|1^]), but we still find it useful to give all the details which we 
are going to use while studying the interacting case. 

The Fourier mode expansion of the V'-operator, appropriate for the boundary conditions 
(H), takes the form: 



]j^^sm{kx)csk (4) 

with k = im/L, n being a positive integer. The single electron spectrum is e{k) and the 
Fermi surface consists of the single point k = kp (see Fig.|I|). Concentrating on the vicinity 
of this point we define slow varying right and left moving fields: 

V2L (5) 



such that (p = rm/L) 

i,s{x) = e*'-^^,B(x) + e-'^'^^^sLix) . (6) 

These fields, however, are not independent, as in the case of periodic boundary conditionslii, 
but satisfy: 

i^siix) = -iPsr{-x) . (7) 

Therefore, one can actually work with the right moving field only, the left moving one is 
then defined by the above relation. The boundary condition 

^.(0) = 

is automatically satisfied, whereas the condition 

ML) = 

implies that the operator ipsnix) should obey 

So we can regard the field ipsnix) as defined for all x but obeying the periodicity condition 
with the period 2L: 

iJsRix + 2L) = xIjsr{x) . (8) 
5 



In terms of the right moving field, the kinetic energy (|T|) takes the form: 

c "J — L 



where we have hnearized the electron spectrum and the energy is accounted for from the 
Fermi energy of a reference system with A'^so number of spin s electrons. 

The right moving Fermi field (^) obeys ^^enoc/ic boundary condition (^), so it can straight- 
forwardly be bosonized. We will conveniently employ the following version of the bosoniza- 
tion formula: 



^-iOs giTT f AAf, ^i4>s (x) 

where ANs is the number of extra electrons with spin s, 

AN, = Ns- iV,o , 
the variable 6,, canonically conjugate to ANg, 

[0s, AN,] = z, 

is defined modulo 2tt. The operators 1], are real (Majorana) fermions, 

{VhVi} = ,r/,^ = 1 , 

which stand to assure the correct anticommutation rules for electron operators with different 
spin s. The phase field 0s (a;) is given by the expression 



(x) = ^ /^e'^'-^^/^^, + H.c. 



q>0 

and satisfies periodic boundary condition: 



qL 



(Psix + 2L) = (Psix) . (10) 

Here bg are canonical Bose operators; g = vm/L, n is an integer, and we have introduced a 
high-energy cut-off a. It is straightforward to check that the operators ipsnix), defined by 
Eq.@, obey standard fermionic commutation relations (in the limit a — 0). 



Alternatively one could write the Fermi operator (^) in a normal ordered form, noticing: 



V2l V2to 
we shall use both normal ordered and not normal ordered forms of Fermi operators. 

The presence of the momentum space cut-off a reflects the finite band-width of the 
original electron band {a is understood to be much larger than the lattice spacing; actually, 
1/a is the region around kp where the electron band spectrum can be linearized). 

The density of right moving electrons is given by 

We notice that 

Psl{-x) = Psr{x) . 
The bosonized form of the kinetic energy is: 

dx : PsRix)psR{x) ■.= VfY. ^blgKg + ——{ANsY . 

sq>0 

Before we turn to the interaction effects, we define the bosonic variables corresponding 
to charge and spin excitations: 

bpia)q = {hq ± , 

and 



B. Interaction effects 

In order to make use of the above bosonization procedure also for the case of interacting 
electrons, we try the same trick; namely, we express the part of the Hamiltonian, which is 
responsible for the interaction, Eq.(|]), in terms of the right moving Fermi field ipsR only. The 
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cost is that the resulting expression is highly non-local in space as illustrated in Fig.^. Nev- 
ertheless, as we show below, at least in the case of short-range electron-electron interaction 
U{x — y), the problem can be quite simply treated in terms of bosonic fields. 

Here we focus on the case of i? ^ a, where R is the characteristic radius of the region 
where the function U{x ~ y) is essentially non-zero. Then the g-dependence of correspond- 
ing Fourier transforms can conveniently be neglected. (This implies that all the distances 
considered should be much larger than R.) 

The interaction part of the Hamiltonian contains several terms classified in what follows. 

The terms, diagonal in the electron densities, can be written in conventional way: 



dx 



Pp{a)R{x)pp{a)R{x) + pf>{^)L{x)pp{a)L{x) 



2 JO 

This renormalizes the sound velocity: 



dxpp(^^)R{x)pp(^^)R{x) . (11) 



, 9p{<t) 

Vp{a) =VF + 



27r 



The term mixing right and left densities is of the form: 



gp{a) dxpp(„)R{x)pp(^„)L{x) = ^gp{a) dxp p{^^)r{x) p p(^^)r{-x) , (12) 

i.e. it is non-local in space (see Fig.||). Still, this term is quadratic in the electron densities 
and therefore takes a simple form in terms of bosonic operators (the remaining terms, which 
do not assume a form quadratic in densities are discussed in the last part of this section). 
Consequently, the Hamiltonian becomes: 



u=p{a) [g>0 



+ 



where 



27r 



This can be diagonalized in a standard way by the Bogolubov rotation 



(13) 



hyq cosh(v9i,)6^g - smh.{^Lp^)hl 



vqi 



(14) 



where 

tanh(2^.) = . (15) 

[Notice that whereas the rotation angles ip,^ are defined in the same way as in the bulk case, 
Eg ■ (p!5D , there is an important difference in sign in Eq.(0).] This rotation is achieved by 
the canonical transformation 

H UHU^ = E hv.qblK + ^(AiV.)4 , (16) 



where 



uq uq ■ AT 



cosh.{2(piy) 

The previously defined v^n can be alternatively written as 



Vi^N = (17) 

being 

= exp{2(p„) . 

The unitary operator U is defined by: 

U = exp\lY. {blqbU - I (18) 

[ ^ g>o J 

The next step is to find how the Fermi operators transform by applying U. Employing 

the method of Lieb and MattisB, after lengthy but straightforward calculations we arrive at 

the main result of this section, i.e. the expression for the electron annihilation operator in 

terms of free bosons for the case of the interacting Fermi system with open boundaries: 

Ui:,R{x, t)U^ = -^e-^'^eir -^^^ (19) 



expji^e 

where e^s is +1 unless s =1 and u = a, when its value is —1. We have defined 
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Cy = cosh(v9,,), Sy = smh.{ipy). 



It is important to notice that, if one would write the expression (0) in the normal 
ordered form, this results in x-dependent pre-exponential factors: 



: exp <^ i ^ 



(20) 



2L 



V2' 



V2' 



where the constant is given by 



and the function P{z) is defined in the Appendix A. Although the pre-exponential factor 
in (^) is x-dependent, it does not depend on t. This explicitates the asjTiimetry between 
space and time coordinates in the present problem (for the translation invariant system, the 
factors resulting from normal ordering are just constant^). 



C. Luttinger liquid picture 

We now pause to discuss how the open boundary conditions modify the standard Lut- 
tinger liquid picture. According to the Haldane analysis^!, a gapless ID system with periodic 
boundary conditions is described by a low energy Hamiltonian of the general form 

H-E'^siblK + ^im' + 'if^r-, (21) 

AA^ being related to the total number of particles, and J to the total current. The different 
velocities in (^) obey the relations: 

vj = Kvs, vn = 

where K is the parameter which governs the asymptotic power low decay of all the correlation 
functions. Two parameters, e.g. vs and vj, are therefore sufficient to determine all the low 
energy properties of the system (concept of the Luttinger liquid universality). For electrons 
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with spin, due to the spin-charge separation, low energy behavior of the system is also 
described by the Hamiltonian of the form ( pT]) in each charge and spin sector. 

As we learned from the above analysis, in an interacting Fermi gas with open boundary 
conditions 

(i) spin-charge separation still occurs; 

(ii) each (spin or charge) sector is described by the Hamiltonian (p!6D, similar to (pi]). 



The important difference with Haldane's analysis is that the Fermi system with open bound- 
aries is bosonized with the help of the right moving Bose fields only, and therefore there is 
only one (in each spin and charge sector) conserved "topological number", AN. The total 
current J is not any more conserved; its dynamics we discuss in Appendix B. The Luttinger 
liquid concept still holds since the two parameters {vs and vn = vs/K) are sufficient to 
describe the low energy behavior of the system. In the spirit of Ref. [|1^], vn and K should 



be treated as phenomenological parameters, though they are related to the interaction con- 
stants g and g which are in turn determined by the interaction potential U{x — y) in the 
starting Hamiltonian (0). 

Notice that we started the above analysis with local (in real space) interactions: gp{x) = 
gp6{x), etc. (The interaction becomes non-local only a/ter formulating the problem in terms 
of right moving fields.) If we would have readily started with long-range interactions of 
radius a (but nevertheless R <^ L), the interaction term in ( pTSD acquires off-diagonal 
corrections (in g-space). In Appendix C we show that these off-diagonal terms do not 
contribute to the asymptotic behavior of the correlation functions for |x| ^ R and VF\t\ ^ R 
and can therefore be neglected. The exponents are thus determined by the g limit of 
the Fourier transform of the interaction constants g^. 

As it is well known, for the case of electrons with spin, the above considered terms - 
quadratic in the electron densities - do not account all the possible interaction processes 
effective at low energy. The remaining process is the so-called spin backscattering process 
which is described by the following term in the Hamiltonian: 
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L 



2 H J_^(^^'^lRi^)^sR{-x)'iplR{-x)'^sji{x) 



(22) 



where s denotes —s. Using Eq.(|T9|) we find that, under the transformation ([T6|), the spin 
backscattering term takes the form: 

rL 

(27ra)2 J-1 

Eq.(p3D can not be diagonahzed in general Nevertheless, one can analyze the effects of the 
spin backscattering by applying the renormalization group (RG) method, as it has been 
done for infinite system0. 

For Kcr > 1 the operator (p3|) is irrelevant, and it fiows to zero under RG process, 
meanwhile renormalizing K„ to a smaller value K* > 1. For spin isotropic interaction the 
fixed point value is 

K: = 1. (24) 

In this case, therefore, the approach we have developed above correctly describes the low 
energy properties of the system^. 

On the contrary, if K^^ < 1, the spin backscattering interaction is relevant, i.e. it flows 
to strong coupling under RG process. This is interpreted as the opening of a gap in the 
spin excitation spectrum. It is known thatS, under scaling procedure, K^,^ decreases and 
eventually crosses the value K^r = 1/2 at which the model has been exactly solved by Luther 
and EmeryEl. Therefore the solution at this point is believed to give a qualitatively correct 
description of the gapped phase for any < 1. The influence of open boundary conditions 
on the gapped phase is analyzed in Appendix D. 

Finally, we notice that although we have modeled the flnite electron system by imposing 
vanishing boundary conditions (^, one could equivalently switch on a binding wall potential. 
For instance a potential V{x) such that V{x oo) but V{x —oo) ^ Vq is equivalent 
(if Vq is larger than the Fermi energy) to a left boundary. In this case one should use as a 
basis the eigenf unctions of this potential ifk{x) while deflning the electron fleld operator: 
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ipsix) = ^^(fk{x)csk- (25) 

Since at large positive x the eigenfunctions ipk{x) are standing waves, the expression (pSj) 
asymptotically coincides with (^) and our above considerations are applicable also to the 
binding wall potential case. 

III. CORRELATION FUNCTIONS 

Given the bosonized form of the electron operators, Eqs.(|T^-(^D|), one can compute 
various correlation functions. Here we discuss some examples of physical interest. 

A. The Friedel oscillation 

Since the system under study is obviously not translation invariant, the mean electron 
density is not homogeneous, so that the Friedel oscillation is build in the ground state. 
Namely, we obtain 

s 

-[P(2x)]^E.^-e-'^(2x) + (x ^ -x) . (26) 



Consider distances x -C L. The Friedel oscillation then takes the form: 



Sp{x) = — - 
7Ta 



a 



sin(2A;i;'a;) 



_Vq^+4x^_ 

This expression describes how the perturbation of the electron density, caused by the bound- 
aries in a semi-infinite system, decays at large distances. 

It is interesting to examine to which extend the open boundary consideration applies to 
doped systems. Kane and Fisher noticed that a local impurity potential 

V{x)J2^lix)^s{x), 
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is equivalent, at low energy, to an infinite barrier, i.e. to a boundary (for infinite system)cl. 
Indeed, they found that the backscattering part of the impurity potential flows to infinity 
under scaling according to: 



where the energy uj is scaled to zero starting from the bandwidth uj = D. Repulsive {Kp < 1), 
spin isotropic {K„ = K* = 1) interaction is assumed. The proof is completed by studying the 
opposite limit of strong potential, which results in a weak link between two independent semi- 
infinite systems - the electron tunneling through a weak link is described by an irrelevant 
operator (which flows to zero under scaling). 

Assume now that many scatterers are present in the system. This is what happens in 
the case of doped quasi- ID electron systems. Does the above argument imply that, upon 
doping, the electron gas will break up, according to the concrete realization of the impurity 
potential, into a set of almost independent segments (to which open boundary description is 
applicable) connected by weak links? To answer this question we notice that these segments 
are finite, so that there is a minimal excitation energy (dimensional quantization), Umin- 
Consequently, the above scaling procedure should be stopped at energies u ~ ujmin- The 
relevant question therefore is whether Umin is so small that the impurity potential grows 
large enough to produce effectively infinite barriers before uo reaches oJmin- Obviously, this 
is just a condition on the scatterers concentration. The energy uJmin should be compared 
with the cross over scale, ujq, between the weak coupling regime of the impurity potential 
and the strong coupling one (the latter being equivalent to a regime of tunneling trough a 
weak link). The latter scale can be estimated by imposing the solution of Eq.(^) to be of 
the order of the bandwidth. That gives: 



Hence, setting Umin ~ vf/xq, where xq is the mean distance between scatterers, we find the 
following condition on doping 




(27) 
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1 /Z}x2/(l-^p) 
xo»7- 77 ■ (28) 



kp Wo, 

Notice also that a short-range character of the impurity scattering was imphed: ~ l/kp, 
where qq is the radius of an individual scattering potential. Otherwise, we should impose 
an additional condition: 

xo > ao . (29) 

Thus, if the conditions (|28|) - low doping - and (|29|) - short-range scattering - are 
satisfied, we arrive at a novel physical picture for doped quasi-lD materials: the impurity 
potential renormalizes to large values thereby breaking up the electron system as shown 
in the Fig.^; the electron tunneling across weak links and from chain to chain can then 
be treated by perturbation theory. This picture essentially differs from those previously 
discussed in the the literature, which deal with perturbation expansions in the disorder 
potentialilli. 

Let us briefly outline some consequences. 



(i) The Friedel oscillation (^) will be frozen in, depending upon the concrete realization 
of doping. Of course, this oscillation will be smeared out by impurity averaging. Still, 
an effect should remain in the density-density correlation function which acquires an 
extra doping dependent 2kp component steming from the impurity averaged Friedel 
oscillation. This should in principle be experimentally measurable (e.g., by neutron 
scattering). 

(ii) The calculation of conductivity is, under these circumstances, an interesting but puz- 
zling task dealing with a type of random resistance network. 

(iii) At low temperatures the electron tunneling between chains becomes relevant, so that 
the 3D character of the problem takes over and a broken symmetry ground state forms 
(typically charge- density wave)0'ili3; the effect of small doping (which increases a phase 
disorder but stabilizes the Friedel oscillation) on the formation of this state is unclear. 
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These issues deserve, in our opinion, further investigations. 



B. Green's function 



The electron Green function is given by (t > 0): 



(30) 



a,b=±l 



where we have substituted Eq.(||) and defined 



Making use of Eq. (|l^) , we find 



GfAx,y;t) 



27ra 



(31) 



where x*^^-* = x ±y, and the phase factor is given by: 



^ix,y,t) = -X 



2 „ 



It is straightforward to calculate the finite temperature version of Eq . (pTD , that is however 
not of an immediate interest. 

Let us now discuss some limiting cases of the expression (|3TD . 



(i) It is interesting to understand how the "bulk" behavior of Green's function is recovered. 
Clearly, one should impose the condition that the relative distance between points x 
and y is much less than all the distances to the boundaries: 



X*- <^ min{a;, y, L — x, L — y} . 



(32) 



However, even if the above condition if fulfilled, the boundaries still influence the 
Green function behavior provided that the time interval t between the creation and 
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the annihilation of the extra electron is large enough to allow the excitations to reach 
one of the boundaries and to be reflected from it. Hence the additional condition: 

{vpt, Vat) <^ min{x, y, L — x, L — y} . (33) 
Provided the conditions (132) and f|33D are satisfied, Green's function flSTF) takes the 



following asymptotic form: 



a 



cl/2 



a 



(~) + Vj^t — ia 



X 



si/2' 



(34) 



which is the electron Green function of the translation invariant system (t > 0). From 
the bulk behavior of Green's function (0) one can deduce, for instance, that 

/1\IE.(^-+^^"') 

{^,{x,t)^l{x,0)) ^ [^jj . (35) 

(ii) On the other hand, if our ^/^-operator is close to the boundary, the condition (^) can 



not be satisfied any more and one should use the general formula (|31|) . Consider the 
case X = y = 0. Of course, strictly at the boundary, a; = 0, the electron operator 
vanishes, Eq.(^). So, by writing x = we mean that x is of the order of a. In this 
case we find: 

GsAO,0;t) Sss' (^jj . (36) 

The latter formula holds actually for any x and y satisfying 

(x, y) < {vpt, v^t) (37) 

and means that the power law exponent of the Green function at the boundary differs 
from the bulk one. This result is implicit in the Kane and Fisher treatment0lll and 
also agrees with conformal field theory consideration^. 

(iii) Different exponents come into play in the static limit G{x,y]t = 0). For t = and 
X = y the expression (^) coincides, of course, with (^^. In the limit x = 0, ?/ ^ a 
(but y <^ L), however, one finds: 
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GiO,y;0)^ i-j - . (38) 

Notice the difference between (|36|) and (^81). 

C. Ward's identity 

Because of tfie novelty of tlie bosonization approach to open boundary systems, the 
question arises whether any important physical processes have been missed in the above 
consideration. Put another way, is it possible to derive the Green function (^Tj), and other 
correlation functions, by means of an alternative method not relying on specific tricks of 
the bosonization approach [like Eq.(^ to Eq.(|l^)]? Alternative method of deriving the 
correlation functions of ID interacting Fermi systems, based on the classification of diagrams 
(Ward's identity), have been devised by Dzyaloshinskii and LarkinEl. Everts and Schulz 
reahzed that Ward's identity can equivalently be derived using equations of motion for 
electron density operatorJ^. 

The starting point, as in Section II, is to formulate the problem entirely in terms of right 
moving fields. The full Hamiltonian then reads: 

= y |wir^/'Jj(x)(-2a^)?/'^(x) + |pij(x)pR(x) + |p/j(x)pij(-x)| . (39) 

(For the sake of clarity we consider the simplest case of spinless fermions.) 
The Dyson equation for the single-particle Green function, 

Gl,{t) = -^{T{c,it)cl,iO)}) , 

where the operator Cp is defined by Eq.(|D with kp + p ^ p, takes the form: 

{tdt - vpp) Gl,{t) = Spp,5{t) -^T. [gF^-.A^, t) + ~gF;\^,{t, t)] , (40) 

being F the two-particle vertex function defined by 

F^At^t') = -^{T{p,{t)cp{t%{0)}) (41) 



with 



P 



The crucial step in solving the ID problem is to find a relation between the vertex 
function F and the Green function G which enables one to express the former in terms of 
the latter thereby resulting in a closed equation for G. Since the non-local interaction in 



(p9|) does not conserve momentum, the diagrammatic approach of Ref . ||19|| is, in our case, 
less convenient than the equations of motion method of Ref. |^ . 
Standard commutation relation! 

[pq,P-q'] = —Sqq' Q, Q' > ^ 

71 

(the 1/2 difference with Ref.0 is due to q taking the values nn/L) lead to the following 
equation of motion for the density operator: 

{idt - vpq) pq = -^qpq + -^qp-q ■ (42) 

Deriving Eq. f^lD in t and substituting Eg . (^21) , we arrive at the version of Ward's identity 
applicable to the present problem: 

{id, - vpq) Fl^,{t, t') = ^qF^Jt, t') + |pgFj,(t, t') ^^^^ 
+i5{t - (t) - t5{t)G^,,,_^{t' - t) . 

The Ward identity ( |43| ) can easily be solved with respect to the vertex function F (in 
cu-space). Substituting this solution into Eq.(^OD and transforming back to the (a;, t)-space, 
we find, after algebraic manipulations, the following closed equation for the Green function 
{t > 0): 

{dt + vpd,) In G\x, y, t) = K{x, y, t) , (44) 

with 

K(x y t) - ~ ^)-^^ _^ (t;^ + v)sc _ vpsc _ (vp - v)c^ _ {vf + v)s'^ 
'' X + y — vt X + y + vt x x — y — vt x — y + vt 
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The limit L ^ oo (semi-infinite system) is taken. Here and below we omit (imaginary) 
short-time cut-off terms (which can easily be restored). Analogously to the case of fermions 
with spin (see Section II), we have defined 



c = cosh(y9), s = sinh((y9); tanh{2ip) 
and the renormalized sound velocity v is given by 



9 



2'kvf + 9 



Vf + 



9_ 

2tx 



9_ 

2tx 



Among many solutions of Eq.(^J), the one which coincides with the free fermions Green 
function for the non-interacting case and is correctly symmetrized with respect to x and y, 
can be written in the form: 



i 


{x + y — vt){x + y + vt) 


sc 


a 


c2 


a 


2'Ka 


Axy 




X — y — vt 




X — y + vt 



(46) 



The Green function (^6]) coincides with the appropriate asymptotic form of spinless ver- 
sion of Eq.(PT|) thus demonstrating the equivalence of "purely" bosonization methods of 
Section II and the approach based on Ward's identity for Fermi systems with open bound- 
aries. [We remind that the total Green function G{x,y,t) is determined by G^{x,y,t) via 
the relation (pO|).] 

The Ward's identity method can straightforwardly be extended to the case of finite L, 
finite temperatures, as well as to the case of fermions with spin. We are not presenting these 
calculations here since they anyway lead to the results identical to those of the bosonization 
technique. 



IV. TUNNELING AT THE BOUNDARIES 

In this section we study several tunneling processes at the boundary of the wire which 
might be relevant for experimental applications. 
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A. Two terminal conductance 



Consider a semi-infinite wire coupled to a normal (3D) metal through an insulating 
barrier, as shown in Fig.^. A potential difference V is applied between the metal and the 
wire and we assume that the potential drop occurs just across the barrier, whereas both the 
metal and the wire have homogeneous chemical potentials. The Hamiltonian describing the 
system can be written in the form: 

H = Hm + Hw+T -VQ (47) 

where Hm{w) is the Hamiltonian of the isolated metal(wire), T describes the tunnehng 
between them: 

^ = H dfT{x,r) \ipwsi^)'>pMsir) + H.c. 

g jw jm '- -' 

where T{x,r) the tunneling matrix element between the wire and the metal, which we take 
limited to the vicinity of the barrier. The last term in (^71) is responsible for the potential 
difference between the wire and the metal, Q being their charge difference 



Q = e{Nw - 



where Nm{w) is the total number of electrons in the metal(wire). 
The potential drop induces a current / which is defined by: 

i(t) = ^ = -2teY: [ dx [ dfT{x,r) kk(^)^Ms(^1 - H.c.l . (48) 

at g Jw Jm -' 

In linear response the average current at time t is given by 

{iit)) = t [dt'e{t-t'){[iit),Qit')])v{t') 



Since [/, Q] = Aie^T, and Q is conserved in the absence of tunneling, the above correlation 
function vanishes if the thermal average is taken over the eigen-states of the Hamiltonian 
(H^) with T = 0. Therefore one needs the next order correction in the tunneling, which 



leads to the following expression 

21 



(/(t)) = z I dt'e{t-t') I dt"e{t'-t"){[i{t)J{t")])oV{t') 



(49) 



where the thermal average (...)o is now taken at T = 0. 

Due to the relation P8|), the current- current correlation function reduces to a product 
of the local single particle Green functions for the metal and the wire. As we have shown 
previously the local wire Green function behaves as 



The latter result is true independently of the presence of electron-electron interaction in the 
3D metal, as follows from the quasi-particle pole residue Z being finite in the Fermi liquid 
theory. 



where Gq is a frequency independent constant which is determined by details of the barrier 
and D is a high energy cut-off related to the metal and wire bandwidths. At zero frequency 
the temperature dependence of the conductance is given by the same formula with the 
temperature replacing the frequency. 

This result enables us to discuss the more realistic geometry of Fig.||. In this case the 
voltage drop V occurs at the two insulating barriers, i.e. V = Vi + V2 (the indices refer to 
each barrier), while the currents flowing across each barrier are the same I = Ii = I2. The 
latter condition holds since no accumulation of charge can occur in the wire. Therefore the 
two terminal conductance 




while the metal Green function is 




Plugging these results in ( ^9[ ) we find for low frequency conductance 
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Thus, in the physically relevant situation of spin isotropic repulsive interaction (/Cp < 1 
and = 1), the two terminal conductance vanishes at low frequency (low temperature) 
with a power law. The conductance of a Luttinger liquid with an impurity has been predicted 
to vanish by Kane and Fisher^. Here we have shown that such an effect should occur even 
in the case of a perfectly clean wire. This reflects the irrelevance of the tunneling T between 
the wire and the normal metal as a consequence of the electron-electron interaction in the 
wire. Thus a two terminal experiment as the one we discussed above can provide a simple 
way of probing the Luttinger liquid behavior in quantum wires. Notice however that those 
wires should be sufficiently long since the dimensional quantization introduces a low energy 
cut-off Vp/L below which the conductance becomes temperature independent and moreover 
gives rise to a finite charging energy (causing Coulomb blockade type phenomena) which we 
have neglected. 



B. Boundary operators 

Similarly to what we have shown for a normal metal - quantum wire - normal metal 
contact, one can imagine other experimental set ups in which different boundary processes 
are involved. In the case we just studied the corresponding boundary operator was the 
single Fermi field operator. For example, at the contact between the quantum wire and a 
standard 3D superconductor the tunneling of electron pairs should be considered^. In Table 
|, worked out by applying Eq. ([T9|) , we give a list of possible boundary operators Ob with 
their scaling dimensions Xb defined by 



V. MAGNETIC IMPURITY 



In this section we consider the effects of a magnetic impurity in an interacting wire. The 
coupling to the conduction electrons is provided by an antiferromagnetic exchange J and a 
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local potential V. On general grounds one expects V ^ J, but even in the unphysical case 
V = 0, a local potential will be generated by the exchange coupling^. The ^ = problem 
has been previously studied by Lee and Toneril and by Furusaki and NagaosaH. In Ref . 



it was found that for repulsive electron-electron interaction the low temperature fixed point 
corresponds to the screened impurity spin, similarly to what happens in the conventional 
single channel Kondo problem for non-interacting electrons. 

Here we study the opposite hmit V ^ J. In this case it is more appropriate to first 
diagonalize the Hamiltonian with only the local potential. For repulsive interaction that 
corresponds to cutting the wire at the impurity site and treating the residual tunneling 
through the barrier as a (irrelevant) perturbationS. Taking into account also the spin ex- 
change, one can in general write the following impurity Hamiltonian: 

i=l,2 a/3 al3 

+ tJ2{4ai^2a + ii-C.), (50) 

a 

where 1(2) corresponds to the right(left) side of the impurity, 5* is the impurity spin-1/2 
operator and a" are the spin-1/2 matrices. The first term represents the exchange interaction 
of the electrons of each lead with the impurity spin. The second one corresponds to tunneling 
processes with spin flip, and the third one to the tunneling without spin flip. 

For r = t = the model is equivalent to the two channel Kondo model, which is known 
to exhibit anomalous behavior. Finite F and t will introduce anisotropy between the two 
channels. For non-interacting electrons this anisotropy is known to be a relevant perturba- 
tion which brings the system back to the single channel behavior at low temperatureS. A 
repulsive interaction {Kp < 1, assuming K^j = 1 to assure spin isotropy) changes the bare 
dimension of F getting it irrelevant, without affecting J (which remains marginal). The 
resulting renormalization group equations for weak interaction are: 



■{r + T 



dT f 1 \ 1 

1 - — F + F J. 



dliaL 271V F 

— = fi- ,-■ 

dlnL \ KpJ irvp 
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Although a small F will firstly decrease under scaling due to Kp < 1, the contemporary 
increase of J will ultimately drive F to larger values. Therefore the relevance of F should 
be more appropriately studied for large values of J. The best way to proceed is then to 
first assume F = 0, let J flow to the strong coupling fixed point, and analyse around it the 
relevance of a small F. 

We bosonize the Fermi fields according to the procedure outlined in section II. We in- 
troduce charge and spin bosonic variables for each lead (piu and (j)2u ('^ = P, cr), and their 
symmetric and antisymmetric combinations (psu and (pau- In terms of these fields the Hamil- 
tonian (EOl) reads 



imp 



2'Ka 
F 



cos 



H.C. 



cos 



H.C. 



ap 



H 5"^ sin(0a^) sin 

TTCt \ JK 



2t 



cos 



COS 



•'ap 



(51) 



In the spirit of the Emery-Kivelson solution to the two-channel Kondo modeS, we can get 
rid of the phase factors involving (psa by performing the canonical transformation 



U 



-iSz4>s 



The transformed Hamiltonian is simply 

1 



imp 



J± COs(0aa) + F_L COS 



■'ap 



sm 



^ap 



+ —S'd, 

271 



2^ X 

H COS(0a^J 

7ra 



cos 



^ap 



(52) 



where X = — 2711) p. If A = F^ = t = the Hamiltonian ( p^ can be easily studied since 
5*^ commutes with the Hamiltonian and can be given a definite value ±1/2. For a fixed 
the Hamiltonian is equivalent to the one describing impurity scattering in a Luttinger liquid 
(for each field (paa and (pap)- If is known that a local cos(/?0)-term is relevant if /5 < -\/2, 
marginal if /5 = a/2, and irrelevant otherwise. Therefore the J± operator is always relevant, 
whereas F_|_ is relevant only if 1/2 < Kp < 1. 
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If both and T± are relevant, the low temperature fixed point behavior is that of the 
single channel Kondo model, as we previously discussed. Just this case has been analysed 
in detail by Furusaki and NagaosaS. Thus, we have shown that the presence of a strong 
potential scattering does not modify their conclusions, provided that Kp > 1/2. 

The new feature which arises from our analysis is that for Kp < 1/2 the two channel 
Kondo behavior is stable with respect to a channel asymmetry. The behavior of the model 
around the two channel fixed point can be performed in a way similar to that of Ref.[ 
p6| . Without repeating the calculation we just remind that the low temperature impurity 
susceptibility x ~ 1^(1/^) the specific heat Cy ~ Tln(l/T). This situation could be 
realized in quantum wires with long-range Coulomb interaction, which may result in a very 



small Kp, see Ref.pT 
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APPENDIX A: 

The computation of average values of exponentials containing phase fields requires the 
knowledge of the sum: 



S(z) = y — e*^^-"'? = In ( — ) + In P{z) + if{z) (Al) 



where : 



P z) = - 



2LWsinh2 (ff ) + sin^ ( 
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fiz) = tan-1 



sm 



Qna/L _ cos 

Let us consider the limit of these functions for z close to one of the boundaries. For \z\ <^ L 
(but \z\ ^ «) we have 



f{z) = tan 



-1 



a 



sign(z). 



The same expression holds at the other boundary {z 2L) with the replacement z — > z—2L. 



APPENDIX B: 

As an interesting application of the open boundary bosonization we consider the following 
problem: let us imagine that in a spinless Luttinger liquid on a ring of length L we insert 
an impurity at the origin. This impurity has two possible states, which can be thought as 
two spin states |t(i))- We model the impurity as a local back-scattering potential^, and we 
assume that the sign of this potential depends on the impurity spin. This amounts to add 
to the electron Hamiltonian a term of the form 

Va,cos{(f)R + (f)L). (Bl) 

where the Pauli matrix acts on the impurity states, and is the phase field corresponding 
to the right (left) moving electrons. Since there is no term in the Hamiltonian which flips 
the impurity spin, az can be given a fixed value, which does not evolve with time. If the 
interaction is repulsive (parametrized by the exponent < 1) we know from Ref. that 
the impurity at low energy effectively acts as an infinite barrier, thus cutting the chain into 
two disconnected leads of length L. Consider now the following x-ray edge type of problem: 
at time r < the impurity is in the state |J,). Suddenly at r = the spin is reversed, thus 
changing sign of the scattering potential, and finally at time t = t the spin is reversed back. 
This process is described by the correlation function: 
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x(t) = ak-(tv+(o)U). 



(B2) 



Can we calculate the long time behavior of (p^ ) ? We recently encountered a similar problem 
analysing the four channel Kondo modefl, and we think it may be relevant also for other 
impurity models. The way to evaluate ( |B2| ) is similar to the standard bosonization approach 
to the x-ray edge singularity. Notice that the unitary operator 



I— J 
2 



U = exp 

where J = Nn — Nl is the total current, has the following property 

U[Ho + V cos(0ij + <Pl)] = Ho-V cos{(I)r + <Pl), 

i.e. it changes the sign of the potential term without modifying the bulk part of the Hamil- 
tonian Hq. Thus U is equivalent to the spin flip operator. This implies that the correlation 



function ( |B2|) can also be written as 



x{t) = (il uit)u\o) li) , 

and in this representation its evaluation is straightforward. Since the low energy fixed point 
corresponds to cutting the ring at the origin, we use the open boundary bosonization to 
rewrite the current operator 

J = Nr - Nl = I dxpnix) - pl{x) = [ dxpn^x) - pr{-x). 
Jo Jo 

By performing the Bogolubov rotation to get rid of the bulk interaction we get 

1 



J 



- 0(0)) 



and therefore 



U = exp 



(0(L) - 0(0)) 



Since the fields are now free Bose fields with logarithmic correlation, we immediately see 
that the dynamics of the total current is characterized by the following correlation function 
{t < L/vf): 
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This result agrees with the analysis of Ref . ||29|| . Hence the correlation function ( |B2| ) takes 
the form: 

X{t) oc (-] . (B3) 

x{t) has the power law behavior typical of the x-ray edge singularity, but the exponent is 
determined by the interaction only and is not bounded from above. 

APPENDIX C: 

This Appendix is intended to show that a long-range electron-electron interaction (i.e. 
interaction of finite radius R) does not modify the conclusions of Section II; namely, that the 
exponents of correlation functions are determined by zero-momentum Fourier components 
of interaction constants. In the case of long-range interaction the bosonic Hamiltonian ([13D 
is no more diagonal in g-space. We proceed in two steps: firts, we demonstrate how the 
diagonal part of the interaction can be separated from the off-diagonal one and, second, we 
observe that the latter is irrelevant (does not contribute to the exponents). 

The interaction term (0) is now of the form 



2 Jo Jo ^^^^^ ~ + Pp(<^)L{x)PpHL{y)\ (CI) 

and the g term (jl^) changes analogously. (The index z/ = p, a is suppressed but the results 
are valid for both charge and spin sectors.) 

Simplifying the problem (but not affecting qualitative results) we work in the limit of 
semi-infinite system [L —>■ oo). Then, in g-space, the Hamiltonian takes the form: 

H = Ho + Hint (C2) 

with the continuum version of the free part. 

Ho = dqvFqb\bq , 
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and the interaction term: 

Hint ■ 



47r Jo Jo ^^^y^ [HQi, (I2) + H-c. 

+ HQi,-Q2) +H.C. 



(C3) 



where 



HQi, Q2) = g{qi, q2) + g{qi, -52) 



and g{qi, ^2) is defined by 



1 /■= 
g{qi,q2) = Re- / 



dp 



(C4) 



27r (p + gi + 2(5) (p - g2 - ^(5) 
Here g{p) is the Fourier integral transform of g{x) and 6 a positive infinitesimal The formula 
( P^ with ^(p) instead of (^(p) defines g{qi,q2)- 

Shifting the integration contour in up in the complex p-plane, we can write: 

i{qu 92) = g{q2)S{q2 + qi) + g{q2)S{q2 - qi) + 52) • (C5) 

The first two terms in (|C5| ) stem from the residual of the pole p = q2 + iS. They are 
responsible for the diagonal part of the interaction. The off-diagonal part, SI{qi, ^2), is due 
to singularities of the functions g{p) and g{p) in the complex plane (see Fig.|^). Since the 
only regular complex function is a constant, the only case when 61 = corresponds to local 
(in real space) interactions. 

In order to illustrate the properties of the off-diagonal part of the interaction we consider 
now a specific example: 



gip) = gip) 



go 



In , 

2 W + A2 



(C6) 



The interaction (CU) may model a screened Coulomb interaction in quantum wires (the 
cut-offs Xg and Aq are then related to the inverse screening length and inverse lattice spacing 
respectively). Evaluating the branch cut integral (Fig.|^) we find: 



Sgiqi,q2) 



go 



E 



(gi + q2) iti 



tan ^ f -T^ I — tan ^ f 



A. 



Aa 



30 



It is important to notice that 51 tends to a constant in the hmit of small q: 

Obviously we can, in general, write: 

^Hqu q2) = goR, for (gi, ga) < l/R ■ (C7) 

Here go is the interaction strength and R is related to the inverse distance to the first 
singularity of g{p) from the real axis (i.e. to the radius of the interaction). 
We can thus re-arrange the total Hamiltonian (|C^) writing it in the form: 

H = H + 5H . 

Here H is the continuum version of the Hamiltonian (|T3|), it contains Hq and the diagonal 
in q part of ( |C3|) . H can immediately be diagonalized [the exponents are determined by 
g{q = 0) and g{q = 0)]. The remaining part, in small-g limit, takes the form: 

goR 



i?2 q2 



The irrelevance of 6H is already clear from its scaling dimension (with is equal to 2). One can 
also demonstrate this more straightforwardly: since SH is quadratic in Bose fields, H can be 
diagonalized. Doing so, we find that the leading corrections (due to 6H) to the correlation 
function {(f){x)(t>{y)) [which behaves as \n{x — y) in the absence of 6H] drop as {R/xY, {R/yY 
at large distances thus do not modifying the asymptotic behavior of fermionic correlators. 

APPENDIX D: 

Here we study how the Luther- Emery solution is modified by open boundary conditions. 
If we define new Fermi operators: 

-\/2vra ' 
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where 9^ is conjugate to ANcr, then exactly at K^r = 1/2 the spin backscattering term 



(E3) assumes a quadratic form in these operators. The Hamiltonian describing the spin 



excitations becomes: 



= V, j ^dx^l{x){-id^)^ ^{x) - i/\ J ^dxsign{x)¥^{x)'^^{-x) , (Dl) 



where 



A = ^ . (D2) 

The sign(a;) factor in (pi|) stems from the phase factor f{2x) in the a — >■ hmit, which is 
necessary to preserve the particle-hole symmetry of the model 

or equivalently 

This symmetry implies that the spectrum should be symmetric around zero energy. 

Being quadratic in the fermion fields, the Hamiltonian ( piD can be brought into diagonal 
form 

where the new Fermi operators are related to the old ones by the canonical transformation 

Cl = I dxXe{x)'^Ux) ■ 

The wave functions Xe{x) satisfy the following non-local Schrodinger equation: 

exe{x) = -iVadxXe{x) - iAsign{x)xe{-x) . (D3) 

Defining the functions x^e^Kx) = Xeix) and xi'Kx) = Xei—x) (for x > 0), we can conve- 
niently rewrite (|D3| ) as the system of differential equations: 
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(D4) 

Those functions obey the following boundary conditions 

Xi^\0) = xi-\0) , (D5) 

Xi-'\L) = xi-\L) , (D6) 

where the latter derives from (|TDp. 

Each particular solution of ( P^ ) has the form exp {±ipx). The dispersion relation turns 
out to be 

e' = Kp)' + A^ (D7) 



and the allowed values of p derive from imposing the boundary conditions ( P5| , P6| ), which 
lead to 

sin(pL)(zA — e — v^p){iA — e + v^p) = . (D8) 

The solutions are p = nn/L with n a positive integer, which correspond to energies |e| > A, 
and p = ±iA/v^, with energy e = 0. Therefore, similarly to what happens in the case of 
an infinite system, the spin backscattering term opens a gap A in the continuum spectrum. 
The novel feature is the existence of the zero energy states, which are bound states localized 
at the boundaries. So, in the L oo limit (semi-infinite system), the wave function of the 
bound state localized at x = takes the form: 

Xoix) = ./5e-AN/.. . 
V v„ 
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FIGURES 

FIG. 1. Single particle spectrum: for the case of open boundaries the Fermi surface consists of 
the single point kp- 

FIG. 2. The system with open boundaries can be described either in terms of right and left 
moving Fermi fields (the upper part of the figure) or in terms of the right moving field only (the 
lower part); in the latter case the density-density interaction becomes non-local. 

FIG. 3. Slightly doped quasi-lD conductor breaks up into independent segments of the electron 
gas (described by the open boundary analysis), these segments are coupled back by weak tunnel- 
ing processes of two types: tperp (tunneling between different chains) and ti2 (tunneling between 
neighbouring segments in the same chain). 

FIG. 4. Contact between a quantum wire and a 3D metal, with an applied potential difference 

V. 

FIG. 5. Geometry for a two terminal conductance measurement. 
FIG. 6. Impurity spin in a quantum wire in the presence of a strong potential scattering. 

FIG. 7. Integration contour in Eq. (|C^) . 
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TABLES 



TABLE I. Boundary operators and their scaling dimensions 



Boundary 


bosonized 


exponent 


operator Oh{t) 


form 


Xb 








V'L(*)V'i?s(t) 


d4>p{t) + d4>a{t) 


1 








'ipRs{t)ll^Rs{t) 
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